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Abstract 



1^ , Extending our recent work [Ij, we report on a calculation of the vector current correlation 

function for light valence quarks in the deconfined phase of quenched QCD in the temperature 
range l.lGTc ^ T < 2.98Tc. After performing a systematic analysis of the influence of cut-off 
^ . effects on light quark meson correlators using clover improved Wilson fermions, we discuss resulting 

CN) ' constraints on the electrical conductivity in a quark gluon plasma. 

o ■ 

(N 1 Introduction 

^ ■ The spectral representation of the correlation functions of the vector current directly relates to the 
>- invariant mass spectrum of dileptons and photons. Additionally in the Umit of small frequencies it 

^ ', determines a transport coefficient, in the case of the vector correlation function of light quarks, the 
^ [ electrical conductivity. 

- - ' At temperatures relevant for current heavy ion experiments non-perturbative techniques are mandatory 
for the determination of these quantities. Perturbative studies of the vector spectral functions [21 13] and 
also the inclusion of non-perturbative aspects through the hard thermal loop resummation scheme |1] 
break down, especially in the low invariant mass region, indicated by an infrared divergent Euclidean 
correlator [5], leading to an infinite electrical conductivity. However it could be demonstrated that 
the spectral function at low invariant masses in fact increases linearly resulting in a finite electrical 
conductivity of the quark gluon plasma ^ . This behavior could also be established by previous work 
using lattice QCD [8J. 

In |T] we analyzed the vector correlation function at T ~ 1.45Tc and performed its extrapolation to the 
continuum limit based Oil precise dcitcL cit various lattice sizes, corresponding to different lattice cutoffs. 
We then obtained reliable results for the determination of the spectral properties and the extraction 
of the dilepton rates and transport coefficients. Here we report on a first extension of this work to a 
wider temperature range and present results on the electrical conductivity in the temperature region 
1.16T, < T < 2.98T,. 
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2 Thermal vector correlation and spectral function 



The desired information is encoded in the large distance region of temporal Euclidean correlation 
functions of the vector current at non-zero temperature and fixed spatial momentum, 

G^y{T,p) = / d^x (J^(r,f)j](0,0)) e'P-"^ , where: J^(r,f) = g(r, f)7^g(r, f) . (1) 



Whereby the current- current correlation functions can be represented in terms of an integral over 
spectral functions, p^yiu^p^T): 

G„(r,AT) = rpM^,p,T) ^"^"'"'[-^^f >' , H = 00,ii,V, (2) 
Jq zvr smh(a;/2T) 

here we denote by pa the sum over the three space-space components of the spectral function and also 
introduce the vector spectral function pv = Poo + Pa- 

The time-time component of the vector correlation function, Gqo{tT), is conserved and can be 
shown to be proportional to the quark number susceptibility, Goq{tT) = —XqT. The vector correlation 
functions Gu and Gy therefore differ only by a constant, 

Gv{tT) = G,,{tT) - x,T . (3) 

At high temperature and for large energies corrections to the free field behavior can be calculated 
perturbatively; the vector spectral function can be deduced from the calculation of one loop corrections 
to the leading order results for the thermal dilepton rate [9J. This yields, 

pvico)^ + co^ tanh(a;/4T) , w/T > 1 . (4) 

27r V vr / 

In the free field, infinite temperature limit the spatial part of the spectral function contains a 6- 
function at the origin. Different from the time-time component, where the 5-function is protected by 
current conservation, this 5-function is smeared out at finite temperature and the low energy part of pa 
is expected to be described by a Breit-Wigner peak [ini [HI IT2] . 

Pii (^) = X9Cw ^2^^j.^^^2 • (5) 

The limit w — )■ is sensitive to transport properties in the thermal medium and in this case it leads to 
the electrical conductivity 

^ = ^lim^. (6) 

T 6 t^^o LoT 

The Breit-Wigner form then yields a{T)/Cem = '^XqCbw / i'^^) ■ Note here that in the infinite tempera- 
ture limit the width of the Breit-Wigner peak vanishes, while at the same time cbw ^, Xg ^ T'^ and 
consequently the electrical conductivity is infinite in the non-interacting case. 

These considerations give rise to a sensible Ansatz describing the high temperature vector spectral 
function: 

Pii{u) = 2x,CBw ^^ + y/2)2 + 1^ + tanh(u;/4T) . (7) 

this Ansatz depends on four temperature dependent parameters; the quark number susceptibility 
Xq(T), the strength {cbw(T)) and width (r(T)) of the Breit-Wigner peak and the parameter k(T) 
that parametrizes deviations from a free spectral function at large energies. 

In the following we will use this Ansatz to extract the electrical conductivity and the vector spectral 
function from the correlator data. 
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2.1 Moments of the vector spectral function 



In addition to the vector correlation function itself we will calculate its curvature at the largest Euclidean 
time separation accessible at non-zero temperature, i.e. at tT = 1/2. The curvature is the second 
thermal moment of the spectral functions at vanishing momentum, 



1 d"GH(rT) 

(jr 



n\ d(rT)' 



1 /"°° do; /tux " ph{uj 



tT=1/2 



(-) ■ , . , H = tt,V, 8 

\l / smh a;/2i 



n\ Jo 27r \TJ sinh(a;/2T) 



where n is chosen to be even as all odd moments vanish. These thermal moments give the Taylor 
expansion coefficients for the correlation function expanded around the mid-point of the Euclidean time 
interval, 



n=0 ^ 



2n 



^-rT] . (9) 



In the infinite temperature, free field limit the integral in Eq. |2] can be evaluated analytically [13] 
and one straight forwardly obtains the first three non-vanishing moments for massless quarks 

^(0),/ree ^ 2^{0),/ree ^ 27-3 Q^Jree ^ ^STT^ ^g ^(4), free ^ 12471^^3 ^^^^ 

3 5 21 

Note here that all thermal moments, G^"^ with n > 0, will be sensitive to the smeared 5-function 
contributing to pa^u) although we expect this contribution to become more and more suppressed in 
higher order moments [1]. 

In the following we will analyze the ratio of Gh{tT) and the free field correlator G^^^{tT)., 

^'^'^^ - ' ^R%^'' - R%^1.^ ("i - ttV + .... 1 , (11 



Gff''''{TT) Q(0)Jree \ ^ \ H HJreeJ ^2 

as well as the ratio of mid-point subtracted correlation functions 

Gh{tT) - Gg) 



A^(rT) 



G'rirT) - G% 



{0)Jree 



^ + (r''''' - Rf^t^ (l-rT\ + .... 1 . (12) 



Q{2)Jree I ^ ^''HJreeJ ^2 

Here we used the notation R^^'™'^ = G^^^ /G^^\ Note that the curvature of these ratios at the mid-point 
determines the deviation of ratios of thermal moments from the corresponding free field values. While 
the ratios of correlation functions differ in the H = ii and H = V channels due to the additional 
constant contributing to Gy(rT), this constant drops out in the subtracted correlation function, i.e., 
Av{tT) = Au{tT) and gJ,"^ = gJ"^ for n > 0. 



3 Temperature dependence of the vector spf on the lattice 

The numerical results presented here are an extension to our earlier study and comprise calculations 
that have been performed at several temperatures 1.16Tc ^ T < 2.98Tc. Consequently all parameters 
were chosen and adjusted following the procedures outlined in [1]. As such they are obtained from 
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Figure 1: The vector correlation function normalized by its free continuum counterpart at varying 
temperature in tT (left) and [fm] (right). 

quenched QCD gauge field configurations generated with the standard SU(3) single plaquette Wilson 
gauge action [H] using the clover improved Wilson action with non-perturbatively chosen clover co- 
efficient csw [IS] in the fermion sector. The temperature range is then varied by keeping the gauge 
coupling fixed for different temporal lattice sizes = 16, 32, 40. Note here throughout we also show 
results at T ~ 1.45Tc for A^.^ = 48 for comparison. The simulation parameters are given in Tab. [H 

In Fig.[T]we show the ratio Ga^rT) / G{['^'^ for the (3 = 7.457 and the (3 = 7.793 (green points) lattices 
in Euclidean temperature units (left) and in physical distance (right). On the right of Fig. [T] we see that 
the (3 = 7.457 lattices lie on top of each other for r < 0.06fm, in view of the results at /3 = 7.793 this 
can be understood as due to cut-off effects. As such the cut-off effects dominate the low distance region 
of the correlation function, in the presented calculations however the cut-off is fixed and therefore we 
expect them to be the same throughout. At r > O.OGfm the deviations should subsequently be due 
to temperature effects. Notice however that the first 6 to 8 points in the correlator are dominated by 
the cut-off and consequently the A^,- = 16 results are more affected by the cut-off effects throughout 
the Euclidean time interval. Nevertheless the = 16 results in Fig.[T] (left) level out in the regime 
tT > 0.4. Note the correlator ratio from the T ~ 1.16Tc lattice exhibits a visibly larger value for all 
Euclidean times, even though its trend is similar to that at T ~ 1.49Tc. 

3.1 Temperature dependence of the thermal moments 

Calculating the temperature dependence of the thermal moments we compute the quantity Ah{tT) 
and then fit to a quartic Ansatz. Using a jackknife procedure we then immediately calculate the ratio 
R^^'^^ / Rf^'^'^'^^'^^ and the corresponding results over T/Tc are shown in Fig.[2l Note here we also give the 
Nj. = 16 and continuum extrapolated results at T ~ 1.45Tc taken from [1] for comparison. Additionally 



Nr 




^ conf 


(3 


a[fm] 


T/T, 


Csw 






Zv 


40 


128 


451 


7.457 


0.015 


1.16 


1.3389 


0.13390 


0.0989(4) 


0.851 


32 


128 


255 


7.457 


0.015 


1.49 


1.3389 


0.13390 


0.0989(4) 


0.851 


16 


128 


340 


7.457 


0.015 


2.98 


1.3389 


0.13390 


0.0989(4) 


0.851 


48 


128 


191 


7.793 


0.010 


1.45 


1.3104 


0.13340 


0.1117(2) 


0.861 



Table 1: Parameters for the calculation of vector correlation functions on lattices of size x Nr. 
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Figure 2: Left: The ratio ^V^i^'"^"^"" 

over T/Tc above and below Tc. Shown are the results for 
/3 = 7.457 (red) and the continuum extrapolation at T ~ 1.45Tc (black). Additionally a number of 
results from A^,- = 16 lattices at varying temperatures (blue) are given, including also results that have 
been obtained using data from [16] 



Right: Table of results for '°V^i 



(2,0),/ree 

a 



here we also show points below Tc from calculations done by the Bielefeld lattice group that have already 
been examined in |16 j. 

Above Tc, i.e. in temperature region evaluated here 1.16Tc < T < 2.98Tc, the data points are seen to 
be almost constant in T/Tc. Turning to the table on the right of Fig.[2]we identify a slight rise of the 
ratio within 0.5% going from ~ 1.49Tc and A^^ = 32 to ~ 2.98Tc and A"^ = 16. 

Below Tc on the other hand we immediately identify a clear temperature dependence of the ratio 
4''°V4''°^'^"'- As such there is a 30% difference between the results at T ~ 0.55Tc and T ~ 0.93Tc. 
At the same time the result at T ~ 0.93Tc is roughly 20% larger than that encountered above T^. 

Together the drop below and the constant behavior above Tc are very interesting from the point of 
view of a possible p-resonance in the spectrum. As such we estimate the impact of a p-resonance on 
the ratio of thermal moments i?!f°V4''°^'^''^' by assuming the p-resonance to contribute a single cosh 
in the correlator: 



G'p(r) = A ■ cosh 



T 



(13) 



taking the second derivative in tT of Gp{T) then gives the second thermal moment, thus we find: 



0(2,0) 



T 



T \ 2 



(14) 



As a consequence we expect the ratio R^p'^^ to go like ~ 1/T^. Adding also a continuous contribution 
to this model ratio we arrive at the following estimate for the temperature dependence below T^. 



d(2,0) _ f 
T<Tc ^cont ~r l 



Cpmp 

Tr 



(15) 



Even though we have only two points below Tc we can nevertheless fit this estimate to the data and 
the result is shown as blue line in Fig.El The qualitatively very different behavior above and below Tc, 
subsequently leads us to conclude that the drop of the ratio R^f^^ /rI"^'^^'"^^^^ below Tc is in fact mostly 
due to a particle contribution. At the same time the constant behavior of the ratio above Tc suggests 
that the p-resonance does not or only very weakly contribute in the examined temperature range. 
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k/Xq 


9/7 n 

X /dof 


l/Cem ■ Cr/T 


2.98 


0.396 


0.93(19) 


3.07(70) 


0.160(19) 


0.54 


0.31(7) 
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0.189(6) 


1.37 


0.30(3) 




0.25 


0.71(2) 


3.75(14) 


0.122(4) 


2.07 


0.24(1) 


1.49 


0.396 


0.98(23) 


3.14(80) 


0.165(23) 


0.01 


0.33(8) 




0.354 


0.97(11) 


3.16(40) 


0.166(11) 


0.07 


0.32(4) 




0.3125 


0.92(4) 


3.47(16) 


0.164(7) 


0.26 


0.31(1) 




0.25 


0.90(2) 


3.53(10) 


0.170(4) 


1.16 


0.30(1) 


1.16 


0.396 


1.00(34) 


3.25(12) 


0.235(38) 


0.01 


0.33(11) 




0.354 


0.99(17) 


3.30(62) 


0.236(18) 


0.01 


0.33(6) 




0.3125 


0.88(6) 


3.89(20) 


0.233(12) 


0.05 


0.29(7) 




0.25 


0.87(4) 


3.89(18) 


0.242(6) 


0.14 


0.29(1) 



Table 2: The fit parameters of tlie Breit-Wigner+continuum Ansatz for T ^ 2.98Tc, 1.49Tc and T ^ 
1.16Tc on lattices with temporal extent A^^ = 16, 32 and 40. Note the fit window is varied between 
TminT = 0.25 and r„j„T = 0.396. . 

3.2 Consequences for the vector spectral function 

To discuss the consequences of our findings on the vector spectral function we invoke the Breit- 
Wigner+continuum Ansatz defined in Eq.[7] and fit to the T ^ 1.16Tc, 1.49Tc and T ~ 2.98Tc data 
respecting also the thermal moments. Additionally we analyze the dependence of the fit parameters 
on TminT, as, with the data being subject to potentially large lattice effects, we expect large errors 
originating from the fit-window we choose. The resulting parameters are summarized in Tab. El 
Here it can be seen that the results at TminT = 0.396 and TmmT = 0.354 have the largest errors, never- 
theless these results should be closest to those in the continuum as only the furthest distance points are 
comparatively free of lattice effects. Without a continuum extrapolation however it is at this point not 
possible to cleanly discern the residual deviation due to the lattice cut-off. Decreasing TminT on the other 
hand leads to visible systematic trends in the parameters, as such the correction term k/xq = kT'^/xq 
and the width T = T/T both increase with decreasing TminT, while the electrical conductivity decreases. 
Additionally the x^ /dof increases with decreasing TminT implying a lower quality fit with lower TminT- 
In light of the argument above these results show that indeed the lattice effects influence the fit more 
strongly with smaller TminT-, as they cannot be taken care of by the employed Ansatz. Consequently 
we attribute the largest part of the systematic trends observed above to the lattice effects. 
Focusing on the parameters themselves note that the correction factor k{T) is the most accurately 
determined. Clearly it shows an increasing trend with decreasing temperature differing by a factor 
~ 1.5 between T ~ 2.98Tc and T ^ 1.16Tc. However as k{T) is coupled to the strong coupling via 
k{T) ~ Os/vr in], this behavior is expected. The same can be seen for the width F and the parameter 
2cBwXq/^ = 2cBwXq/ ^T), they also slightly increase with decreasing temperature. However here, 
they remain within errors of each other and the increase between T ~ 2.98Tc and T ~ 1.16Tc is small. 
Naturally this also leads to only small deviations of the electrical conductivity. Note however that the 
results at Nj. = 16 differ by ^ 5% from those of = 32, while this deviation is only ~ 2% for A^,- = 40. 
With the Nt- = 16 lattice being most affected by the systematics of the fit, we assume the largest part 
of this deviation to find its origins there. Even so we find the electrical conductivity to be to a good 
degree constant and within ~ 7% across the temperature region evaluated here, whereby the values at 
the individual temperatures have fit errors around ~ (7 — 11)%. 
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In the temperature range 1.16Tc ^ T < 2.98Tc Our results therefore imply that the electrical conduc- 
tivity depends linearly on temperature as expected from perturbation theory p] and at /3 = 7.457 we 
find the value: 

a= (0.33 ±0.04 ±0.11) -Cg^-T , (16) 
whereby the first error quoted corresponds to the temperature effects and the second to the fit errors. 



4 Conclusions 

Extending our earlier study we presented results on the vector correlation function in the high tem- 
perature phase of quenched QCD in the temperature range 1.16Tc < T < 2.98Tc. First results show a 
similar behavior to those extensively analyzed at T ~ 1.45Tc and even at T ~ 1.1 6Tc we find no sizable 
contribution attributable to the p-resonance. 

Even though these calculations could not yet be extended to include also a continuum extrapolation as 
in [T], we were nevertheless able to extract the relevant spectral properties to calculate the electrical 
conductivity via an Ansatz and found this phenomenologically relevant transport coefficient to depend 
linearly on temperature with a value of a = (0.33 ± 0.04 ± 0.11) ■ Cem ■ T. 

Given the exact relation between the electrical conductivity and the diffusion constant, a = Xg -D, it 
would be interesting to compare these results with those obtained in the case of charm diffusion [17] or 
from the momentum diffusion constant k calculable via HQET [T8|[T9] and we will do so in an upcoming 
pubhcation. 
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